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Abstract 

We set up and analyse a model for the non-equilibrium evolution of a network 
of cosmic strings initially containing only loops and no infinite strings. Due to 
this particular initial condition, our analytical approach differs significantly from 
existing ones. We describe the average properties of the network in terms of the 
distribution function n(£, t)d£, the average number of loops per unit volume with 
physical length between I and I + d£ at time t. The dynamical processes which 
change the length of loops are then estimated and an equation, which we call the 
'rate equation', is derived for dn(£,t)/dt. In a non-expanding universe, the loops 
should reach the equilibrium distribution predicted by string statistical mechanics. 
Analysis of the rate equation gives results consistent with this. We then study 
the rate equation in an expanding universe and suggest that three different final 
states are possible for the evolving loop network, each of which may well be realised 
for some initial conditions. If the initial energy density in loops in the radiation 
era is low, then the loops rapidly disappear. For large initial energy densities, 
we expect that either infinite strings are formed or that the loops tend towards a 
scaling solution in the radiation era and then rapidly disappear in the matter era. 
Such a scenario may be relevant given recent work highlighting the problems with 
structure formation from the standard cosmic string scenario p|. 
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1 Introduction 



Cosmic strings are just one of a number of different topological defects which may have 
formed at a phase transition very early in the history of the universe @, || |J. Since 
they are topologically stable, cosmic strings may have survived until today and thus may 
provide an alternative theory to inflation for explaining the formation of structures such 
as galaxies in the universe, and the temperature fluctuations in the cosmic microwave 
background (CMBR). Over the last twenty years, a 'standard' picture has emerged of the 
initial distribution of strings produced at the phase transition: numerical simulations of 
string formation || suggest that this distribution consists of two distinct 'populations'. 
The first are cosmic string loops (strings which are closed on themselves) whose length 
distribution is smooth, typically a scale invariant £~ 5 ^ 2 power-law. The second is a sep- 
arate population of infinite strings (we define the difference between a long loop and 
an infinite string below). Numerical simulations || [7|, |[ and analytic models J9|, [T(J of 
the subsequent evolution of such a network in a critical expanding universe then suggest 
that the strings reach a 'scaling' solution in which the energy density in strings is a fixed 
fraction of the energy density of the universe and all scales in the problem grow as the 
horizon. Given this 'standard' picture, calculations of the observational effects of cosmic 
strings together with comparison with experimental data suggest that GUT scale strings 
may not be responsible for both the fluctuations in the CMBR as well as the density per- 
turbations [0] . These authors suggest that strings or other active perturbations might be 
reconciled with data if they could somehow disappear rapidly at a redshift of about 100. 
Alternatively one could consider the effects of strings in a universe with a cosmological 
constant |TJ, |TJ. 

However, one should recall that many very important details are still unknown both 
about the formation and evolution of cosmic strings. For example, there is uncertainty 
in the order of early universe phase transitions JT3) : can topological defects form at all in 



continuous transitions? If they can, what would their initial distribution be and would 
it still contain the two different populations? Even in first and second order transitions, 
analytic verification of the numerical results of string formation is very difficult given the 



non-equilibrium and non-perturbative nature of the problem [14| . Recently attention has 
once again turned to understanding the effects of the lattice |15], [16], [T7|, 0, [L^, ^(| |21 |. 
There are also obvious limitations to numerical simulations of the evolution of cosmic 
string networks. For instance, because of the complexity of the system, these simulations 
can only run for a small fraction (~ 10~ 25 ) of the total time between string formation and 
today. Also important gravitational back reaction effects have yet to be incorporated in 
numerical work and any string structure which is smaller than the lattice size cannot be 
resolved. These limitations mean, for example, that the small scale structure on strings 
is not treated as accurately as it needs to be. Analytical models for the evolution of 
a network of cosmic strings containing infinite strings have also been developed. One 
such model has been proposed by two of us || , and in it we included the effects of back 
reaction. The production of loops from infinite strings was also treated, but we ignored 
the subsequent fate of the loops assuming that they all decayed by gravitational radiation. 

In this paper we develop a rather different analytical model for the evolution of a net- 
work of cosmic strings that contains only loops. Before outlining some of our motivation 
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for considering this problem, we define what we mean by infinite strings and loops. The 
reason for doing this is that there might be some confusion about the difference between 
very long loops (whose length tends to infinity) and infinite strings. Also, most definitions 
of infinite strings are tailor made to strings on lattices with periodic boundary conditions 
B. Here we work instead in an infinite space — an Q = 1 universe — and hence there 
are no boundary conditions. To define the meaning of an infinite string, consider a cube 
of side X, and let n(£) be the number density of closed strings with length £ which we 
assume to be Brownian random walks of step £. (We want to show that by definition n(£) 
must only describe a loop distribution without infinite strings.) Since the length density 
L = J °° £n(£)d£ in strings must be finite, then if n(£) oc £~ a for large £ this implies that 
a > 2. Within the cube, the total length in strings which do not touch the walls is 

L m = J \X-yJin £n(£)d£. 
By definition, the total length of closed strings in the volume is given by 

roo 

Uot = X 3 / n(£)£d£. 
Jo 

If, in the limit X — > oo, the ratio Lj„/L tot — > 1 then we define this to mean that there 
are only loops and no infinite strings; that is, as X —>■ oo no strings are long enough to 
cross the boundary. Since this ratio is indeed 1 for all n(£) for which L is finite, it follows 
that the distribution n(£) does describe solely loops. Infinite strings take the form of a 
separate distribution which is characterised by an X independent length density L in f] 
this just adds a contribution to L tot so that typically the ratio L in /L tot — > 0.25 [|, |22"fl . 
For any X, working with only a loop distribution means that there are less strings which 
cross the box than there would be if a population of infinite strings were also present. 

Our motivation for considering the evolution of such a network of loops comes from 
different directions. Firstly, it is interesting to focus on the loop network alone (this, we 
believe, has never been done), and hence try to understand its properties better. For 
example, it is often said that infinite strings are responsible for seeding the perturbations 
which then collapse to give galaxies. However, could a loop network not survive long 
enough to perform the same task? Existing numerical work suggests that the evolution 
of a loop network may differ significantly from that of one that also contains infinite 
strings |23] . Loop networks may therefore have very different cosmological consequences. 



In particular, might the evolution be such as to 'save' cosmic strings — i.e. that the loops 
scale in the radiation era and then disappear in the matter era QXJ] ? Finally, as a result 
of a very different dynamical simulation of a first order phase transition, Borrill [16[] has 



argued that there is no evidence for infinite strings, but only a population of loops. As 
yet this claim has not been disproved even though it has been contested |T7|, [Tj| |2"T| . It 
is certainly the case that the fraction of infinite strings decreases if there is a variation in 



domain size at the phase transition |L9| . Given our previous comments on the complicated 
nature of phase transitions, we do not believe that an initial distribution containing only 
loops is ruled out for certain. 

In this paper, we therefore start from the hypothesis that the initial distribution of 
cosmic strings formed at the phase transition contains only loops and no infinite strings. 
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We set up a model for the evolution of such a network and analyse some of its properties. 
Our model can be applied both to a loop network evolving in a non-expanding universe, as 
well as to a network in a critical Q — 1 expanding universe. In the non-expanding case, our 
dynamical approach gives results consistent with those of string statistical mechanics. In 
an expanding universe there are no results to compare with, and our goal is to determine 
whether the loops ever reach a stable scaling solution. 

This paper is organised in the following way. Section |2] introduces the model and 
contains a description of our approach. We also specify how loops interact with each 
other and derive an equation for the rate of change with time of n(£,t)d£. This equation 
is subsequently referred to as the 'rate equation'. In section [3] we estimate some of the 
coefficients which appear in it and specify our approximations and assumptions. In an 
expanding universe, it is useful to change variables to 'scaling' variables and this is done 
in section f|. Sections and |6] are concerned with the large £ solutions of the rate equation, 
both in a non-expanding as well as an expanding universe. Questions of stability and 
scaling are also discussed. We summarise our conclusions in section |7|. 

2 The Model 

In this section we set up a model to describe the non-equilibrium evolution of a network 
of cosmic string loops. Such evolution may result from the expansion of the universe, or 
alternatively it could describe the evolution of the network towards an equilibrium state 
in a non-expanding universe. Throughout we use a statistical description of the loop 
network. 

2.1 Boltzmann equation for a loop network 

The network of relativistic cosmic string loops is described by the distribution function 

n(£,t)d£, (2.1) 

so that n(£,t)d£ is the average number of loops per unit volume with physical length 
between £ and £+d£ at time t. (For an analogy with a system of particles, see [^4], p6||). 
The dimensions of n are thus given byQ 

Note that the loop distribution n(£, t) is taken to be homogeneous and isotropic in space 
so that it does not depend on the position r of the loop.0 

An equation for dn(£, t) / dt can be obtained once length-changing interactions between 
the loops have been specified. In order to do this, we assume throughout this work that 
when two different loops intersect, or when a loop self- intersects, partners are exchanged 
||. !|]. Ignoring the effects of gravitational radiation and expansion for the moment (see 
below), the two dominant length-changing interactions between loops are then 

1 We set c = 1. 

2 The distribution function could depend on other parameters such as the mean square velocity of the 
loops. Here we consider the simplest case and ignore all dependence other than that on I. 
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A loop of length £ can intersect another loop of length £' to give a loop of length 
£ + £'. 



• A loop of length i can self-intersect to produce two daughter loops of length £' and 
£-£'. 

Consider the first of these processes. If the initial lengths of the two loops are in the 
range I — > I + d£ and £'—>£' + d£', then the number of collisions per unit time per unit 
volume of this type is 

A(£ + £';£,£')(n(£,t)d£)(n(£',t)d£'). (2.3) 

As expected, this is proportional to the number of loops with each of the initial lengths 
times the probability A(£ + £'; £, £') of the collision occurring. Equation ( [2.3|) defines A 
which will then be estimated in section |3|. Dimensional analysis gives 

[A{£ + £';£,£')} = (Time) 2 . (2.4) 

Observe that by definition A must be symmetric; 

A{£ + £';£,£') = A[£ + £';£',£), (2.5) 

and that this process can both remove loops of length £ and produce them (as a result 
of the intersection of two smaller loops). 

Similarly, the number of daughter loops of length £' produced from the self-intersection 
of a loop of length £ — > £ + d£ per unit time per unit volume is given by 

B{£ -£',£'■ £){n{£,t)d£). (2.6) 

The function B{£ — £', £'; £) defined by ( |2.6|) has dimensions 

[B(£ -£',£';£)] = j^- 2 . (2.7) 

(lime)" 1 

Again note that B is symmetric with 

B{£ -£',£';£) = B{£', £-£';£). (2.8) 

We will estimate both A and B in section |^. 

In a complete microscopical description of the system, these A and B processes would 
be related though time reversal. Here, however, we are considering the average properties 
of the network and label the loops only by their length. Since this does not constitute a 
complete set of state labels for the loops, we no longer expect A and B to be related in 
any simple way. 

2.2 Derivation of the rate equation 

An equation for dn/dt can be obtained in terms of the functions A and B. Recall first 
that in an expanding universe, £ is not a constant between collisions, £ ^ 0, since the 
length of the loops changes for example though gravitational radiation ^ [| . Let H be the 
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expansion rate, and consider a length interval from £ to £ + di. In a time At, the change 
A(n(£, t)d£) in the number of loops per unit volume with length £ — > £+d£, is given by the 
number of loops per unit volume which scatter into the interval minus the loops per unit 
volume which scatter out. Taking care with the special interval from to di, the resulting 
rate equation is obtained by letting At — > so that A(n(£,t)d£)/ At — > (dn/dt)d£. It is 



dn 
~dt 



,„ ,d£ dn{£,t) • v 

-n £,t — )^ 1 £+ hm 

di d£ e^o+ 



+ 
+ 



1/2 



A(i; £', £ - £')n(£', t)n(£ - £' , t)d£' 
B{£,£! -£; £')n{£',t)d£' 



- n(£,t) / A{i + i';i,i')n{i',t)di' 



-3Hn(i,t) 


(2.9) 


Term 1 


(2.10) 


Term 2 


(2.11) 


Term 3 


(2.12) 


Term 4. 


(2.13) 



n(£,t) t' 2 B( £',£-£' -£)d£' 
Jo 



In line (|2.9| ), the first two terms can be written as the single differential —d{n£)/d£. Hence 
they represent the net positive flux of loops in £ space. The third term in ( |2.9| ) guarantees 
that 

n(£,t) = W<0, Vt, (2.14) 

as is required for a physical distribution of loops. (For a discussion of this point see also 
section |5.1| .) The final term in ( |2.9| ) arises from the expansion of the universe. 

The four scattering integrals ( |2.1(J| )-( |2.13 ) are subsequently referred to as 'term 1' to 
'term 4' respectively. (Terms 3 and 4 are defined without the minus signs.) Term 1 gives 
the rate of production of loops of length £ from the intersection of two smaller loops of 
length £' and £ — £'. Notice that because of the symmetry of the A coefficient ( |2.5| ), the 
upper limit in term 1 is £/2 so as to avoid over-counting. Term 1 gives a non-linear and 
non-local contribution to the rate equation. Term 3 is the converse of term 1 and is the 
rate at which loops of length £ disappear due to intersections with other loops. Loops of 
length £ may be produced as a result of a larger loop self-intersecting — this is term 2. 
Its counterpart is term 4 which gives the rate at which loops of length £ disappear due 
to self-intersection. Again the upper limit of term 4 is 1/2 so as to avoid over-counting. 
These four terms together with the first part of the equation mean that the rate equation 
is a non-linear and non-local integro partial-differential equation: hence difficult to solve. 
In fact, it is interesting to note that in existing discussions of loops (see for example ||]), 
the effects of terms 1 to 4 are never included — that is, the interactions between loops are 
ignored (though this is not a very realistic approximation). The resulting rate equation 



is then easy to solve, as we shall see in section |5 

In a non-expanding universe we set £ = H = so that dn/dt is just given by the four 
scattering terms ( |2.10|) - (|2.13| ). The rate equation therefore loses its partial differential 
aspect. An important check on this equation is to verify that the total length density L 
in loops defined by 

roc 

L = I n(£,t)£d£ (2.15) 
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is preserved.^ That is, in a non-expanding universe one should have 

dL r°° .dn(£,t) , 



dt Jo dt 



0. (2.16) 



By definition each of the A and B scattering processes individually preserve length (sec- 
tion |2.1| ). However, starting from the rate equation ( |2.13| ) in a non-expanding universe, 
it is in fact difficult to prove (|2.16|) as one often has to deal with diverging integrals. For 
example, if the distribution n(£,t) ~ £ -5 / 2 for large £, then the contribution to ( |2.16D 



from each of terms 1 to 4 is divergent. These need to be carefully regularised in order to 
proceed, but once this is done a long calculation shows that ( |2.16| ) is indeed satisfied as 
required p7| . 



We intend to solve for dn/dt. First, however, the functions A, B and £ must be 
estimated. 



3 Determination of A, B and £ 

We begin by summarising our assumptions about the length scales expected to charac- 
terise the loops. 

3.1 Length Scales 

Throughout this work we assume that there are two length scales on the loops (see 0). 
The first is the persistence length, £(t), the distance over which the loops are correlated 
in direction. Both in non-expanding as well as expanding universes, we assume that the 
loops are Brownian random walks on scales much greater than £. The second scale is 
denoted by ((t) and it measures the 'wiggliness' of the string due to the kinks produced 
as a result of intercommutation fl, |9H . 

As yet we have not been able to obtain equations for t;(t) and ((t). The following 
assumptions are therefore made about their behaviour: 

• In a non- expanding universe the two scales are taken to be time independent; £ = £ c 
and C = Cc where £ c and ( c are constants. 

• In an expanding universe they are assumed to 'scale': 

C(*) = Cc*. (3.1) 

Typically, we expect £ to be of the order of the horizon size whilst the smaller scale 
C should be characterised by the scale of gravitational back reaction. Hence for 
GUT strings we expect ( ^ TG/i^ ~ 1CT 5 ||. 

Expressions for A, B and i may now be found in terms of £, ( and I. Our aim will 
then be to try to solve ( 2.13) for n(£, t) and determine whether, in an expanding universe, 



the energy density in loops also scales. In a non-expanding universe we should obtain 
results consistent with string statistical mechanics. 

3 Length density L and energy density E are related by E = [iL, where \i is the energy per unit length 
of the string. In the following we often use the two words interchangeably. 
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3.2 Coefficient A 



By definition, A(l\ + £ 2 ; £1, £2) is th e probability that a loop of length £1 intersects a loop 
of length £ 2 to form a loop of length t\ + £ 2 (see equation (|2.3| )). We therefore expect A 
to take the form 

A(4 + 4;^i,4)=X^ 2 , (3.2) 

where x is a relative velocity times a geometric factor — see 0. 

For a motivation of this expression, consider two cylinders of lengths T\ and r 2 much 
greater than their diameter d (so that terms of order djr\ and <i/r 2 can be neglected). 
Then, by averaging over all possible orientations between the cylinders, one can show 
that their scattering cross section is proportional to T\Ti times their relative velocity. If 
we model the loops as being made up of £/£ negligibly thin cylinders of length £, then 
the scattering cross section for two loops of lengths £\ and £ 2 must be of the form given 
in (|3.2|) . Note also that ( |3.2|) has the correct dimensions from ( |2.4j) . As in |§, we expect 
X ~ 0.1 - 0.2. 



3.3 Coefficient B 

This coefficient is much harder to estimate. However, consider a loop of length £ and ask 
for the probability per unit time per unit volume of it decaying to give a loop of length 
y with y < £/2. (From (|2.8| ), we then know B for all y.) 

Firstly suppose that £ <C y <C £/2. On these scales the loop is, by assumption, a 
Brownian random walk, and we take B to be given by 



where \ is another relative velocity. To obtain equation ( |3.3|) , consider one of the £ 
segments of the loop, and another such segment m = y/£ steps away. Then the straight 
line distance between these two segments is m 1//2 £ (assuming Brownian random walks), 
and the volume contained between them is of order (m 1//2 £) 3 . Thus the probability that 
these two segments intersect in time dt to form a loop of length y = m£ is equal to the 
ratio of the volume swept out by one of them divided by the total volume: 

^ dt (3.4) 



( m l/2f)3" 

Here x is the relative velocity of the two segments. Since there are £/£ pairs of such 
segments, one obtains (|3.3|) . Again as in [[J, we might expect x ~ 0.1 — 0.2 but it is not 
necessarily equal to x as it is a different relative velocity. 

The case for which y <C £, £ is rather more difficult to estimate. Numerical simulations 
H suggest that intercommutations lead to a rapid build up of kinks on these small 
scales, and therefore B{£ — y, y; £) might be expected to peak about y ~ £. (However, if 
particle production is the dominant mechanism by which the loops radiate energy — as 
has been claimed in |^0[ [28| — then there will no longer be any small scale structure. This 
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interesting case is not considered here.) The B coefficient should also be proportional to 
the length £/£ of the loop and have dimensions given by ( |2.7|) . We therefore write 



B{l-y,y)l) = ^ 2 (j\h( y ,C,£) (</«£)• 



(3.5) 



Here h yy,(,£,) is a yet unspecified dimensionless function which is probably peaked about 
y ~ ( and satisfies h(0, (, £) = 0. 

Given these forms of A and B, we can use the principle of detailed balance and observe 
that if an equilibrium situation is reached in which the rate at which loops break into two 
daughter loops, £ + £' £, £', is equal to the rate at which they join via the A process, 
£ + £', then the corresponding distribution n (£) must satisfy 

A{1 + £'; £, £')n (£)n (£') = B{£, £'; £ + £')n {£ + £'). (3.6) 



For large £ equations ( |3.3| ) and ( |3T2|) show that ( |3.6|) is satisfied by distributions n (£) oc 

e -/3^-5/2_ 

In this paper we focus on determining the large £ asymptotic behaviour of solutions to 
the rate equation, and it appears that the trends reported in sections |5] - |] are in fact not 
sensitive to the form of h (see |[27|). For this reason, we simplify calculations by working 
with the following function B(y, £ — y;£) which we assume to be valid form all £, y: 

B(v£ v . £) - m-yW + v?' 2 (37) 

Here rj is an undetermined length. First note that by definition B is symmetric as required 
( ^78|) and that it has the correct dimensions (fX^) . Secondly, for £ ^> y ^> r/ (|3.7|) correctly 
reduces to ( |3.3| ), whereas for y ~ 7] £ it takes the form ( |3.5| ) with h = y^ 3 ^ 2 (y + ri)~ 5 / 2 . 
Given that ( p.7|) is also peaked about y = £ — y ~ r], we might expect r] ~ (, but we 
shall leave rj arbitrary. However, as with the other lengths characterising the loops, we 
assume that rj scales in an expanding universe with rj = at (where a is a constant which 
should not be confused with the a coefficient introduced in section |2.1| ) and rj is a constant 
in a non-expanding universe. Observe that, given ( |3.7| ), the corresponding equilibrium 
distribution n (£) is 

3.4 Estimating t 

Expansion of the universe, redshift and gravitational radiation (G.R) change the length 
of a loop. We do not consider the loops in the friction dominated era of the universe 
though if we did, an extra term would have to be included to I P, 3, fl0[|. Also, the 
possibility that loops may decay predominantly by particle production is not taken into 
account |2(J j28fl . 

The contribution to £ from redshift and expansion is H 



expansion 
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H (l - 2(v 2 ) e , c ) £. (3.9) 



Here (v 2 )i^ is the average velocity squared on the loop. As a first approximation, we take 
(v 2 )i£ to be a constant with 1 — 2(v 2 )i^ = a. From [|7| we might expect that a ~ 0.1 in 
the radiation era, and a ~ 0.26 in the matter erasfj The contribution from gravitational 
radiation is modelled by 

Here T is a constant, and / (|) is a function which tends to 1 as £ tends to zero, and for 
large £ is proportional to £/(. The reason for this latter limit is that the dominant source 



of gravitational radiation from long loops is from the kinks |[29|| , and the ratio £/( just 
counts the number of kinks on a loop of length £. The former limit comes from applying 
the quadrupole formula to a loop in Minkowski space [§. The simplest form of /, the 
one we assume, is / (J?) = 1 + m| where m is a constant. Of course it will be important 
to check the sensitivity of solutions of the rate equation to this choice of /, and also to 
the approximation made for B. 
Summarising, we write 

K£ 

£ = -TGfi + — , (3.11) 

t 

where 

K = -^ + ^. (3.12) 

Cc 2 

Given the values of a, we expect K < 0.1. For long (many times the horizon size) loops, 
we also expect the dominant length changing effect to be from the expansion of the 
universe, and therefore that K > 0. Another reason for taking K > is given in the next 
section. 

In a non-expanding universe, we set TGfi = K = so that £ = 0. 

Given these estimates for the coefficients A, B and £, the rate equation may now be 
analysed. However, before doing this a final change of variables is carried out so as to 
clarify the problem in an expanding universe. 



4 Rate equation in scaling variables 

In the previous section we have made the assumption that in an expanding universe, the 
two scales £ and ( are proportional to t. Since we aim to see whether the energy density 
in loops also scales in that case (i.e. L oc t" 2 ), it is useful in an expanding universe^ to 
change variables from £ and n(£,t) to 'scaling' variables, z and g(z,t) defined by 

, = { (4.1) 

and 

g(z,t)=t 4 n(£,t). (4.2) 



4 However, since there are no infinite strings here, we expect a actually to be rather smaller than these 
values which will then represent upper bounds. 

5 This change of variables does not make sense in a non-expanding universe. There one must continue 



to work with (2.13) in which I = H = 0. 
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Notice that z is dimensionless, and for lengths which 'scale' it is a constant. Also observe 
that 

fit 2 

where E sc is the energy density in scaling variables. Thus the energy density in loops 
scales when g(z,t) = g(z); these are the scaling solutions we are searching for. Now, 
using the definition of the coefficients A and B, one has that 



poo 1 fOO f] (f\ 

Lit) = / £n(£, t)di = - / zg(z, t)dz =: -=^2 (4.3) 
Jo t z Jo ut z 



A(zt;z't,(z - z')t) = t 2 A(z;z',z - z'), (4.4) 

1 
f 2 



B(zt,{z' -z)t-z't-U) = -B(z,z' -z;z';Q, (4.5) 



where for the moment we have exhibited £ as a fourth argument. Thus a rate equation 

for g(z,t) can be obtained by substituting into ( f2.13j ) and using equation ( |3.11| ). It is 

t^M = g{Zj t) ( 4 -2-*)+2?M (*(!-*) + TG») (4.6) 

- TGfi6(z)[g(z,t)]\ z ^ (4.7) 

rz/2 

+ / A(z;z',z- z')g(z',t)g(z- z',t)dz' (4.8) 
Jo 

/•oo 

+ / B(z,z' - z;z')g(z',t)dz' (4.9) 



- g(z,t) / A(z + z';z,z')g(z',t)dz' (4.10) 

Jo 

- g(z,t) f" 2 B(z',z- z']z)dz'. (4.11) 

Jo 

Here the parameter p is defined by the expansion rate; R oc t 1 ' p , so that p = 2 in the 
radiation era, and p = | in the matter era. Observe that the RHS of the above equation 
is not explicitly t dependent and hence that it is consistent to have g{z,t) = g(z) - i.e. a 
scaling solution is consistent with the equation. It remains to find the form of this scaling 
solution and to see whether it is stable. Henceforth we always work with the scaling 
variables z and g{z,t) in an expanding universe. 

An equation for dE sc /dt can be obtained by using equations ( |4.3|) and ( |2.16| ), and 
integrating (^4.11|) . It is 

dE sc ( T ^ 3\ E sc TGu 2 N sc 

— —= [2 + K-- — " — -. 4.12 

dt \ p) t t v ' 

Here N sc (t) = t 3 N(t) where N(t) = J °° n(£,t)d£ is the total number of loops at time t. 
Thus it is possible for the energy density to scale both in the radiation and matter eras 
provided K > 0, and TGfi ^ 0. We therefore take K > 0, whilst the important role 
played by TGfi will again become apparent in section [|. Note that the assumption of 
scaling lengths £ and £ does not automatically imply that the energy density in loops 
also scales in an expanding universe: this will only happen if values of E sc and N sc are 
obtained such that the right hand side of ( |4.12j ) vanishes. The aim of the rest of this 
paper is to see whether such scaling solutions exist and whether they are stable. 
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The work presented here is entirely analytical and given the complicated nature of the 
rate equation we do not attempt to find any exact solutions but instead focus on their 
large z ^> £ asymptotic behaviour. We are currently trying to develop methods to solve 
(Ell|) numerically |0]. 



5 Search for solutions g(z, t) = z a e 



-f3z 



In the rest of this paper we focus on the functional form of large z solutions of the rate 
equation. This section in particular tries to determine whether stable scaling solutions 
of the form 

g(z,t) = c(t)e- mz z a{t) 0>O, (5.1) 

exist for large z>^ c . 

The motivation for looking for such solutions comes from results of string statistical 
mechanics |3~T| , |32| , |33fl : in a non-expanding universe the loop network evolves until it 
reaches a stable equilibrium distribution, and if the energy density in loops is less than a 
critical value, this distribution is given byf] 

ce~^ z 

fc) = ^sjr ( 5 -2) 



for large z BT], . This form cannot of course be correct for z — > since it would lead 



to infinite energy. Hence it must be modified for small z. In fl5.2|) (3 > since the energy 
in loops must be finite. At (3 — the energy density in loops is equal to the critical value 
and the loops have a scale invariant distribution g(z,t) oc z~ 5 ^ 2 . A phase transition then 
occurs, and the system is unstable against the formation of an infinite string |31, 32, 33]. 



Since our model of loop evolution can be applied to loops in a non-expanding universe, 
it should also predict stable solutions of the form fl5.2|) as long as the energy density in 
loops is below the critical value (which in turn should be determined by the parameters 
in the model). In this section we therefore use these comments as motivation to look for 
large z solutions of the rate equation of the form (|5.1|) both in a non-expanding as well as 
an expanding universe. Here we only consider the case for which /3 > 0. Critical power 
law loop distributions with (3 = are analysed in the next section. 

This section contains a number of parts. In the first we try to get a feel for the rate 
equation by considering a simplified, but unphysical, limit of this equation. In this case 
we show that the equation admits solutions of the form ( |5.1| ) for large z. The conditions 
required for these solutions to be scaling solutions, g(z,t) = g(z), are also specified, and 
comments are made about stability. The aim of section |5.2| is to determine whether a 
distribution of the form ( |5.1| ) is consistent with the full rate equation. To do that, the 
RHS of the rate equation is calculated for distributions (|5.1D, and consistency conditions 
allow equations for a, (3 and c to be obtained. These are then analysed in sections |5.3| - 



5^4] . In section [O] we work in a non-expanding universe. We see that as long as the 
energy density in loops is below a critical value, E crit — which is expressed in terms of 
the parameters of the model — there is a stable scaling solution of the form (|5.2| ). The 



corresponding value of c is determined, and we show that if the energy density is greater 



6 Recall that in a non-expanding universe z = I and g(z, t) = n(£, t). 
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that E crit , there is no stable solution (|5.2|). The results of statistical mechanics suggest 
that infinite strings may be formed in that case. An expanding universe is considered in 
section [|]4] and we investigate whether a scaling solution of the form (|5.1|) can exist in 
the radiation and matter eras. 



5.1 Conditions for a solution g(z,t) = c(t)e f ] ( t ) z z a ( t ) 

We begin the search for large z solutions of the rate equation ( [4.11| ) by considering a 
simple, unphysical, limit of that equation in which there are no interactions between the 
loops; A = B = 0. In the radiation era p = 2 the rate equation ( |4. 11| ) now reduces to 

9g(z,t) /5 \ dg(z,t) dg(z,t) 

t — — — =g(z,t) ( g — J U ~ K ) z dz + TG V—fa TGii5(z)[g(z 1 t)]\ z ^o+- 

(5.3) 

Without solving (|5.3|) explicitly, one can investigate whether it has solutions of the form 
( |5.1|) for large z and all times by simply substituting the ansatz g(z,t) = c(t)e _/3( -') 2 ,2 a O 



into ( |5.3j ). For z > it gives 

t- - tz/3 + talnz = (- - K + a(l - K) - TGfip) - zf3(l - K) + ^E^. (5.4) 
c V2 / z 

In the large z limit where \TGfia(t)/z\ is negligible compared to all the other terms, one 
can equate the coefficients of z°, z and In z in Q5.4Q to obtain equations for a, (3 and c. 
They are 



a = 0, (5.5) 
t$ = {1-K)P ^/3(t) = /5(t )(f) 1 K , (5.6) 



t- = --K+{l-K)a-/3TGfi. (5.7) 



c 



Generally in this section we require /3 > 0. However, in this special case we may note that 
if P(to) = then (3 stays zero for all times. In that case a scaling solution (d = /3 = c = 0) 
can only exist for the special value a = —(5/2—K) / (1—K). Such power law distributions 
are discussed further in section ||. 

If P(to) > then (3 grows with time (as 1 — K > 0) and so the system never reaches 
a scaling solution. From equation ( |5.7| ) observe that when (3 ~ 1/TG/i, the dominant 
term on the RHS of that equation is given by (3TGn , and so c — * as t —>■ 00. Thus 
g(z,t) — > as t — >• 00 and if this were a realistic model of loop evolution, we would say 
that all the loops disappear. The solution (3 = is therefore an unstable one. 

This analysis is confirmed by calculating the exact solution to ( |5.3|) . For an initial 
distribution at t = to given by g(z,t ) = c(t )e~ ,3oz z~ 5 ^ 2 , it is 

g(z,t) = c(t ) ( - ) ^ -^d(z). (5.8) 
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We comment that the effect of the term involving a 5 function in ( [5.3| ) is the appearance 
of the 9 function in the solution ( |5.8|) . Hence the S function ensures that g(z,t) = 
\/z < 0,Vt as was claimed in section [2.1| : had the 5 function not been there, then there 
would be no 9 function in ( |5.8|) and hence the solution g(z, t) would have been unphysical, 
since for any t > t , g(z,t) ^ when —TGfi(t/to — 1) < z < 0. 

We stress that as long as (3 ^ 0, the factor TGfj, plays a very important role in 
determining the large t behaviour of c(t). This is consistent with the comments of section 

§ 

Finally it is interesting to note that if there was an additional term Pg(z,t) on the 
RHS of ( |5.3p that did not vanish for j3 — 0, then [3 = would no longer be a solution to 
the equation. For negative P the solution would instead be an unstable one with (3 ^ 0. 



5.2 The full rate equation 

Now consider the effect of adding the four scattering terms to the RHS of ( |5.3| ). However, 



rather than working with a single value of a as in ( |5.1|) , we work with a more general 
distribution function 

(z + afl 2 ^ (z + a) n 



This can be regarded as a modification of ( |3.8|) . We wish to see if such distributions are 
consistent with the rate equation, and if so to obtain equations for the c n and (3. Note 
that if there is a solution g(z,t) = c§e~^ z {z + a)~ 5 / 2 then we expect the c n — > (n > 1) 
as t — > oo whilst Co(t) tends to a constant. 

Equations for (3 and c n are obtained by substituting the ansatz (|5.9|) into the full rate 
equation, expanding in powers of -, and then equating the coefficients of z, z°, z" 1 . . .. 
As an example of how to deal with the scattering integrals, consider term 2 which is given 

by 



Term 2 = e~ Px V c n / y -— r^dy. (5.10) 

We are interested in the z 3> a limit as we are searching for large z solutions to the 
rate equation. If (3 > (the case considered in this section), the exponential inside 
the integrand naturally provides a cutoff at x ~ 1//?, so that provided we look at large 
z 3> 1//?, the denominator of (|5.10 ) can be Taylor expanded in powers of |. The result 



is 



oo 



Term 2 = ^ c -* £ " 

£c ^ n=0 



1 1 
+ «- + «^ + ... 

z Z A 



z^a,-j, (5.11) 



where each • is an n dependent coefficient arising from the Taylor expansion. The case 
for which (3 = and n ^ (and hence such a Taylor expansion is not permissible) is 
considered in section ||. 

After a similar analysis for the other scattering integrals, equations can be obtained 
for $ and c n . To zeroth order, we set c n = for n > 1 and obtain equations for (3 and c 
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from the coefficients z and z° respectively in the rate equation (ill ). They areQ 

4 



t/3 



c 



(l-K)f3 + ( X c -x) 

/3 3 3/ 
V 2 " P + ~2 



h 



3a 1 / 2 



3 3 3K \ 

o - - + -S- - r ^/3 - 8v^(xc - x), 



where 



ii 



ye 



o (y + a) 5 / 2 



d y = e Pa Jpr{-l/2,/3a) - ae f3a p 3/2 T{-3/2, pa). 



(5.12) 
(5.13) 

(5.14) 



To first order (where c n = 0, n > 2) we have a more complicated set of equations valid 
for P > 0; 



*/3 

t6\ 



>l-K)P + { X c -x) 



h 



4 



3a 1 / 2 



+ XCi 



4 



15a 3 / 2 



( 3 3 , 3i H a s 1/2/- ^ i 8 X c oCi 

V 2 ~ p ~2~ j °° ~ TG ^ C ° ~ 8c ° a (*Cq - X) + 5-^172"' 

+ ^) Cl - rG/i + /? Cl ) + ^a 3 / 2 c (xco - X) 



(5.15) 
(5.16) 



1 3 5K 
2~ p 

3C! 



+XCi 



28c a 1/2 ) + X(8a 1/2 ci - C1J3) + £(1 - K)ac . (5.17) 



Here 

and 
h-- 



3a 1 / 2 
ye" 



,1/2 



„ to+0 )r/.*'= e *( 1+? r ^ 



r(-|.^«)-5^75. (5-18) 



ye 



(y + a) 5 / 2 



dy = e pa \p- 1/2 T(l/2, Pa) - 2ap 1 ' 2 T(-l/2, Pa) + a 2 /5 3/2 r(-3/2, pc 



(5.19) 



5.3 Analysis in a non-expanding universe 



Consider first the behaviour of equations ( |5.12|) -( |5.17|) in a non-expanding universe where 



we should obtain results consistent with string statistical mechanics. For the purposes of 
an analytical analysis, we assume that P is small. Then the zeroth order equations (|5.12 ) 
and ( |5.13| ) reduce to 



P = -2V^^P( X co-x), (5.20) 
c'o = -8y/a(xc - x)c - (5.21) 

Observe that a stationary solution with c = $ = exists when c = Cq™ 4 where 

<%* = X/X- (5-22) 

7 We set £ c = 1 in the rest of this paper so as not to have to write explicitly factors of £ c . These can 
be reinstated into the following equations by changing x — ► 
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As expected, the value of (3 is not fixed by the value of Co in this case: these two variables 
are linked by the total energy density of the system which is conserved by the rate 
equation (section |2T2j ). Also observe that c > x/x =^ c 'o < 0, $ < and similarly 
c o < x/x =^ c 'o > 0, $ > so that this solution is a stable one - as is required for 
equilibrium distributions. 

Figure [l] shows trajectories of constant energy density in the (cq,/3) plane. We 




crit 



crit 
crit 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 



Figure 1: The (cq,(3) phase diagram showing lines of constant energy (solid curves). The 
parameters are a = 1 and x — X — 1 so that Cq™* = 1. The line of stable solutions is 
given by the dashed curve. 

observe that as long as the energy density E is less than a critical value E CT u given by 
(here c = c C Q lt and (3 = 0) 

E *« = f l jdwr* dz = lira (5 - 23) 

then the stable scaling solution is reached since the trajectories of constant energy cross 
the stable equilibrium line for which Co = c c q %1 . If E>E crit then (3 — > 0, but at (3 = 
the equations ( |5.20| ) and ( |5.21| ) break down (see section |]). However, we recall that from 



string statistical mechanics [^T[] the point (3 = corresponds to a phase transition and 
the formation of an infinite string. If this occurs, we cannot analyse the evolution of the 
system further as our equations only model a network of loops which does not contain 
infinite strings. This case is considered in section |J 

Our conclusions on the evolution of (3 and Co remain unchanged by the analysis of the 
more complicated equations ( p.l5| )-( ^7T7D containing c , C\ and (3. Indeed, we see that 

16 



Ci — > for all initial values of the parameters whilst /? and Co behave as discussed above 
(see figures § and ||). 



. 999999 
log t 0.999998 
0. 999997 
0.999996 



0.025 
-0.05 



Figure 2: Evolution of (3, Co and C\ according to equations (|5.15|) - (|5.17|) adapted for a 
non-expanding universe. The parameters are \ = x — 1 so that Cg™* = 1, and a = 1 
giving a critical energy of E cr a = 1.33. The initial values of the variables at time t = 1 
are (3(1) = 0.05, c (l) = 0.9, ci(l) = 0.1 giving a constant energy of E — 0.887 < E crit . 
Since the initial value of Co is less than Cq™*, both Co and /3 are seen to increase according 
to figure (|5) while ci — > 0. 




Figure 3: This figure is similar to the previous one except that now (3(1) = 0.01, c (l) = 2, 
Ci(l) = 0.1. Hence E = 2.57 > £? cr ^. As expected, (3 decreases until it reaches zero, c 
also decreases whilst c\ eventually tends to 0. At (3 = the equations are no longer valid. 

To conclude, in a non-expanding universe we obtain results consistent with those of 
statistical mechanics provided the total energy in the loop network is less than E crit . 

5.4 Analysis in an expanding universe 

In an expanding universe, the small (3 limit of equations ( |5.12|) and ( |5.13|) is 



t(3 = (1 - K)(3 - 2^(3( X co - x), (5.24) 
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t c A = (* _ I + 3 K \ _ - 8v^(xc - x). (5.25) 

c V 2 P 2 / 

The scaling solution is given by $ — c'o = 0. In the radiation era these two conditions are 
compatible — we will discuss the stability properties (or lack thereof) of this solution in 
the next subsection. In the matter era, however, no full scaling solution is permitted. The 
reason for this difference between the two eras can be seen from ( |5.24|) -( |5.25| ). Imposing 
the conditions c'o = $ = gives 




(5.26) 




If TGfi/3 is very small compared to the other terms (we expect TGfi ~ 10 5 and (3 to be 
initially small), the condition for a scaling solution to exist is therefore 

> 0. (5.27) 

Thus in the radiation era (p = 2), a scaling solution can exist as K > 0. In the matter 
era, however, condition ( |5.27| ) requires K > 1/3 > 0.33. Since we expect < if < 0.1, 
this therefore implies that no scaling solution is possible in the matter era. 
Equations ( |5.24| )-( p.25| ) are now analysed in the radiation and matter eras. 

5.4.1 Radiation era 

In the radiation era, scaling occurs when (3 = j3 cr u and Co = c^™* where 

IK 

<**"-» = TWS' (5 ' 28 > 

9ttK 2 

Pcnt ~ Ua{l-Kf (5 ' 29) 



Here we have ignored terms of order TGfi/3. In this limit, equation (|5.25 ) can be written 



as 

t- =8xv^(cS"*-c )- (5.30) 

Co 

Hence Co = Cp™* is an attractor since Co > c^ lt =^> c'o < and Co < Cq %1 =^ c'o > 0. 
However, analysis of equations ( |5.24j ) and ( |5.25| ) shows that, depending on the initial 
values of (3 and Co at time t = tj, there are two very different possible scenarios for the 
evolution. 

Consider first the case in which c (t/) = c§ %t . If /3(tj) > (3 C rit then the dominant 
term in (|5.24j) is the first one, and so (3 increases. Even though equations (|5.24j) and 
( p.25| ) are only valid for small (3, the full equations ( |5.12[) and ( |5.13| ) show that (3 keeps 



on and on growing. When (3 ~ 1/TGfi, then the dominant term on the RHS of ( |5.25|) 
is —YGn(3 and so for large times, c — > 0, (3 — > oo: the loops disappear. If (3 < (3 crit 
then we see from (|5.24j) that (3 — > whilst c — > c^ lt . The scaling solution (3 = (3 crit 
and Co = Cg™* is therefore an unstable one. For arbitrary values of co(£/) and /3(tj), 
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Figure 4: The (cq, (3) plane at the initial time t = tj in the radiation era. The parameters 
are x = X — 1> an d a = 1 so that the critical point, marked with a dot in the figure, is 
at Cq" 4 = 1.0187 and f3 crit = 0.00545. A scaling solution with (3 = j3 crit and c = Cq" 4 is 
reached for initial values of co(£j) and P(ti) lying along the curve. However, this curve 
is unstable, so that for initial values deviating slightly from these /3 — > oo; Co — > or 
P — > 0; Co — > Cq 7 * 4 * as indicated. The effect of decreasing a is to increase the region of 
parameter space for which (3 — > 0. 

the situation is summarised in figure ^. Comparison of figures [| and |^ shows that in 
both cases, (3 —>■ for high loop densities. However, for low densities the loops reach 
an equilibrium configuration in a non-expanding universe but disappear in a radiation 
dominated universe. 

Analysis of equations ( |5.15| )-( [5TT7j ) shows that the effect of adding the c\ contribution 
is to change slightly the position of the line of scaling solutions in figure § but other 
than that not to alter the behaviour of (3 and cq. Furthermore in all parts of the plane 
c\ — ► (see figures | and||). Finally, we observe that when (3 — > oo it grows with time as 
(3 ~ t 1 ~ K . Correspondingly c and C\ tend to zero exponentially; c , C\ ~ exp (— !-fr|^ — )• 
On the other hand when (3 — > it decreases as Pit) 1 / 2 ~ P{ti) 1 ^ 2 — Aln(t/tj) where 
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Figure 5: Evolution of (3, c and C\ in the radiation era according to equations ( |5.15| )- 
( gl7| )- Here (3(1) = 0.05, c (l) = 1.2, ci(l) = 0.01, rG/i = 10~ 3 and the other parameters 
are as in the previous plots. 




Figure 6: Evolution of j3, Co and c\ in the radiation era with (3(1) = 0.01, Co(l) = 1.2, 
Ci(l) = 0.01. The other parameters are as in the previous plots. 

5.4.2 Matter era 



In the matter era, we have argued that there is no scaling solution. Equation ( |5.25|) 
though can still be written in the form (|5.30 ), where now 



8v^(xcr - X) = - Q - ^) <0 (5.31) 

since K < 0.1. In this case therefore Cq"* may be positive or negative. We see that c^ nt > 
if and only if y/a > ^~ (h — ¥f)- As in the radiation era, Co — > c^ lt . On the other hand, 
if eg"* < then c -> 0. 

Let us first consider the case in which Co(t/) = Cq"* > 0. Then since (xCq"* — X) 
is negative, ( p.24|) shows that (3 is always positive and so (3 grows. It increases until 



(3 ~ l/TGfjL at which stage Co — > 0. We therefore see that in contrast to the radiation 
era, the case for which co(tj) = c^ %t always leads to disappearing loops. 

However, independently of the value of Cq™*, it is possible to find values of co(ij) 
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for which (3 — > 0: from ( |5.24j ) this occurs when Co(t/) is very large (see figure |7|). The 
corresponding region of parameter space though is much smaller than in the radiation 
era, so the loops are much more likely to disappear in the matter era. It seems that 
rather extreme initial conditions are required to ensure /3 — > 0. The overall situation is 
summarised in figures [7|-§. In all cases, analysis of equations fl5.15j )- (|5.17| ) shows that 
ci -> 0. 




Figure 7: The {cq, (3) plane at the initial time t = tj. The solid line is the same curve as 
in figure [| for the radiation era and a = 1. The dotted points are related to the matter 
era (again a = 1). Here, in contrast to the radiation era, there is no possible scaling 
solution for all times: though Co tends to a constant as long as a 1 / 2 > — ^-J, there 

is no solution with f3 = as t — > oo. Therefore the points on the dotted line do not 
correspond to (unstable) scaling solutions as they did in the radiation era. Instead this 
curve just marks the boundary between the regions for which (3 — ► (on and above the 
dotted line), and (3 — > oo. In all cases addition of C\ shows that c\ — > always also in the 
matter era. 



6 Search for power-law solutions 

We have seen that both in an expanding as well as in a non-expanding universe, if the 
initial value of Cq is large enough the solution will evolve towards (3 = 0. What happens 
when this point is reached? 
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Figure 8: Evolution of (3, c and c\ in the matter era according to equations ( |5.15| )- 
( |5.171 ). In all of these figures co(tr) = 1.1 and a = 1. From figure 0, the value of (3{ti) 
which cuts off between (3 — > and (3 — > oo is /3 C "* = 0.0000725. In the first figure 
= 0.0004 > /3 c "*(t/) and /? -> oo. In the second /?(*/) = /3 crit (t 7 ), whilst in the third 
P(tj) = 0.00007 < /5 crii (t/). 



In section 572, we performed a Taylor series expansion in inverse powers oia/z, looking 



for a solution of the form (|5.9| ). However, this expansion is invalid when (3 = 0. There the 
integral ( 5.10| ) representing Term 2 is no longer an analytic function in the neighbourhood 



of z = oo. It can be expressed in terms of hypergeometric functions, and expanded in 
powers of (a/z) to give 



Term2 = ^ E ° 



1 1 

Z 1 ' 2 z 



[z > a). (6.1) 



The appearance of half-integral powers implies that the ansatz ( |5.9| ) fails at this point. 

Physically, this failure is unsurprising. We know that in equilibrium in a non-expanding 
universe, there is a phase transition, the Hagedorn transition, at the point (3 = 
f32| , |31| , [2~2]| . There is a critical value of the energy density, E cr i t . If we start with a 
low density, E < E crit , the system will tend towards an equilibrium state described by 
fixed values of (3 and c , where the curves cross in figure [I]. For large density, however, 
the equilibrium state contains not only a loop distribution with g(z) oc z-^ 2 (i.e., (3 = 0) 
but also a single infinite string. The normalisation of the loop distribution (i.e. Co) is 
independent of the density; as the density increases more and more of the length goes 
into infinite string. 

Our equations as they stand do not explicitly allow for infinite strings. The transition 
is signalled by the fact that the equations break down at that point. If we start in a 
non-expanding universe, with E > E crit , the solution will evolve toward (3 = 0. It will 
reach that point in a finite time. However the resulting value of cq exceeds the maximum 
value for which a stable solution is possible, c c ^ lt . What we believe happens therefore is 
that the excess energy is transferred into the form of long strings, leaving a scale invariant 
loop distribution characterised by c = d^ %t . 

At first sight it may perhaps seem surprising that a pure loop distribution can evolve 
to create infinite strings. However, it is not difficult to see that this can indeed happen. 
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Figure 9: The evolution of c corresponding to the previous figure. Thus c (ij) = 1.1, 
and we have the different values of (3(ti) as given above. In the first two figures cq will 
eventually tend to zero when (3 gets very large. 

Creating an infinite string obviously requires an infinite number of intercommuting events, 
but the number per unit length or unit volume is finite. If we start with two infinite 
sinusoidal strings approaching one another transversely, they can break up into an infinite 
set of loops (see figure Conversely, an infinite set of loops can collide to form a pair 
of infinite strings. Indeed if we start with a sufficiently large density of loops, we must 
expect that mutual collisions will lead to a percolating system, i.e., to infinite strings. 

It should be possible to augment our equations by allowing also for a certain density 
of infinite string. We have not attempted to do this in detail (because of certain technical 
difficulties). However, qualitative examination of the problem suggests that infinite string 
could not be formed or survive if the loop distribution corresponded to (3 > 0; in such 
a case the probability of creating loops from infinite strings would be infinitely greater 
than that of creating infinite strings from loops. For (3 = the equations would appear 
to be consistent only when Cq = Cg™*. Thus infinite strings must always be accompanied 
by the invariant loop distribution known to be present in equilibrium. 

We believe essentially the same results would apply in a radiation or matter dominated 
universe. If the initial string distribution consists only of loops, but with a density above 
some critical value, they will evolve towards a distribution in which infinite strings are 
present ((3 = 0). However, such an initial distribution is perhaps rather unlikely: if the 
density is large enough, then in all probability the initial distribution itself would already 
contain infinite strings. Physically, the more interesting situation is the one in which 
the initial density is below the critical value, so that the distribution contains only loops 
for all times. It will then evolve as described in the previous sections and the loops will 
disappear in the matter era. 

7 Conclusions 

We have set up a model for the non-equilibrium evolution of a network of relativistic 
cosmic string containing only loops, taking into account the possibility that loops can 
join to form longer loops, and that a loop can self- intersect. The effects of gravitational 
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Figure 10: Two sinusoidal infinite strings travel towards each other generating an infinite 
set of loops. 

radiation, redshift and expansion on the length of a loop were also considered. 

In a non-expanding universe we showed that the model conserves the total energy 
density in loops, E oc J °° n(£, t)£d£, as required for consistency. Furthermore, under the 
assumption that the two length scales £ and ( characterising the loops are constant in 
time, we were able to show that if E is less than some critical energy density E cr i t , the 
loops reach a stable equilibrium distribution of the form n(£) = CQ Ttt £~ 5 ^ 2 e~ /3e . The value 
of E crit is determined by (3 = and Cq"* = x/X't that is by the ratio of two dynamical 
parameters in the model. For E > E crit we argued that f3 — » and hence that a scale 
invariant distribution of loops is formed. Furthermore, by combining the fact that energy 
is conserved with an extension of the model to include infinite strings, we argued that 
when the point f3 = is reached a distribution of loops of the form n(£, t) = 
formed with the rest of the initial energy density going into infinite strings. These results 
are consistent with those of statistical mechanics. 

In an expanding universe we showed that the situation for the evolving loop network 
is very different. In a radiation dominated universe a scaling solution was observed to 
exist, but it was seen to be an unstable one. Also, depending on the initial energy density 
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in the loop network (figure we saw that either (3 — > oo so that the loops disappear, or 
(3 — > 0. In a matter dominated universe, there is no possible scaling solution. Given the 
analysis of section § we conclude that in an expanding universe, the network may evolve 
in the following way — depending on the initial distribution of loops formed at the phase 
transition: 

• If the energy density in loops is small, the loops disappear both in the radiation 
and the matter eras. 

• If the energy density in loops is greater than a critical value then (3 — > in the 
radiation era, and we may expect the following possibilities: 

1. The point j3 = is reached still in the radiation era or conceivably in the 
matter era, in which case infinite strings are formed with a scale invariant 
distribution of loops. The rate equation then breaks down and the effect of 
infinite strings need to be included (see for example ||). 

2. If (3 does not reach zero we expect the loops to disappear in the matter era. 

We thus believe that the only way in which the loop network can evolve so as to stay as 
a loop network (that is one not containing infinite strings) is for the loops to disappear 
in the matter era. This scenario may be interesting from the view point of structure 
formation from cosmic strings 

Throughout this work we have assumed that in an expanding universe £ and ( scale. 
A full analysis would not make this assumption but itself derive equations for £(£) and 
C(t), though this would be very difficult (c.f. 0). Also, we have analysed the evolution of 
specific initial distributions g(z,t) oc z a e~ f3z and not arbitrary ones which are outside the 
scope of our analytical work. However, given that results from string statistical mechanics 
and simulations of cosmic string formation predict initial loop distributions of the type we 
have analysed, we believe we have studied the most relevant case. The study of arbitrary 



distributions requires a numerical approach [30] 



There are a number of directions in which we would like to extend this work. Firstly we 
hope to study carefully the effect of the transition from radiation to matter domination, 
and hence to clarify which of the above outcomes is the most likely. We also hope to 
explore the cosmological consequences of the evolution of such a network and calculate 
for example the gravitational radiation emitted. If infinite strings are proved to exist 
then we hope to add their effect to this model, for example by combining it with the 
3-scale model |§. 
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